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ISSUE:  Numerical  models  for  predicting  the 
initial  fate  of  material  disposed  in  open  water 
are  required  for  the  following  activities: 

•  Address  environmental  concerns  related  to 
the  disposal  of  dredged  material. 

•  Provide  input  for  long-term  sediment 
transport  models  used  in  disposal  site 
management. 

RESEARCH:  Dredged  material  disposal 
models  were  developed  under  the  Dredged 
Material  Research  Program  (DMRP),  1973- 
1978.  Under  the  Dredging  Research  Program 
(DRP),  additional  developments  to  the  earlier 
models  have  resulted  in  the  numerical  dis¬ 
posal  model  called  STFATE  (Short-Term 
FATE)  for  application  to  split-hull  barge  and 
hopper  dredge  disposal  operations.  These  de¬ 
velopments  have  been  guided  by  both  field 
data  and  data  from  large-scale  laboratory 
tests.  In  addition  to  guiding  model  develop¬ 
ments,  data  from  the  laboratory  tests  have 
been  used  in  model  validation  efforts. 

SUMMARY:  A  numerical  model  called 
STFATE  (Short-Term  FATE)  for  computing 


water  column  concentrations  and  bottom  depo¬ 
sition  resulting  from  a  dredged  material  dis¬ 
posal  operation  has  been  developed  as  a  result 
of  extensive  modifications  to  an  earlier  dis¬ 
posal  model.  For  example,  the  bottom  col¬ 
lapse  phase  that  occurs  for  most  disposals  has 
been  redeveloped  based  upon  an  energy  con¬ 
cept.  These  developments  have  been  guided 
by  large-scale  laboratory  tests  of  disposal  op¬ 
erations  from  a  split-hull  barge  and  a  multiple 
bin  disposal  vessel.  Results  from  the  disposal 
tests  have  also  been  employed  in  model  vali¬ 
dation  efforts.  Comparisons  of  niodel  results 
with  these  data  are  presented  in  the  report. 

AVAILABILITY  OF  REPORT:  The  report 
is  available  through  the  Interlibrary  Loan  Ser¬ 
vice  from  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  Library,  telephone 
number  (601)  634-2355.  National  Technical 
Information  Service  (NTIS)  numbers  may  be 
requested  from  WES  Librarians. 

To  purchase  a  copy  of  this  report,  call  NTIS  at 
(703)  487-4780. 
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Preface 


This  study  was  authorized  as  part  of  the  Dredging  Research  Program  (DRP) 
of  Headquarters,  U.S.  Army  Corps  of  Engineers  (HQUSACE),  and  was 
performed  under  the  "Numerical  Simulation  Techniques  for  Evaluation  of 
Short-Term  Fate  and  Stabilization  of  Dredged  Material  Disposal  in  Open 
Waters"  Work  Unit  32465.  This  work  unit  is  part  of  DRP  Technical  Area  1 
(TAl),  Analysis  of  Dredged  Material  Placed  in  Open  Water.  Research  was 
conducted  during  the  period  September  1988  to  September  1994. 

Messrs.  Robert  Campbell  and  Glenn  R.  Drummond  were  DRP  Chief  and  TAl 
Technical  Monitors  from  HQUSACE,  respectively.  Mr.  E.  Clark  McNair,  Jr., 
Coastal  Engineering  Research  Center  (CERC),  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  was  DRP  Program  Manager  (PM)  and  Dr.  Lyndell 
Z.  Hales,  CERC,  was  Assistant  PM.  Dr.  Nicholas  C.  Kraus,  Senior  Scientist, 
CERC,  and  Dr.  Billy  H.  Johnson,  Hydraulics  Laboratory,  were  Technical 
Managers  for  DRP  TAl.  Dr.  Billy  H.  Johnson,  was  the  Principal  Investigator 
for  Work  Unit  32465. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K  Howard,  EN. 

Additional  information  can  be  obtained  from  Mr.  E.  Clark  McNair,  Jr., 
Program  Manager,  at  (601)  634-2070  or  Dr.  Billy  H.  Johnson,  Principal 
Investigator,  at  (601)  634-3425. 


The  contents  of  this  report  are  not  to  he  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

cubic  yards 

0.7645549 

cubic  meters 

feet 

0.3048 

meters 

inches 

2.54 

centimeters 

square  feet 

0.09290304 

square  meters 

Summary 


A  numerical  model  called  STFATE  (Short-Term  FATE)  for  computing 
water  column  concentrations  and  tx)ttom  deposition  resulting  from  a  dredged 
material  disposal  operation  has  been  developed  as  a  result  of  extensive  modifi¬ 
cations  to  an  earlier  disposal  model.  For  example,  the  bottom  collapse  phase 
that  occurs  for  most  disposals  has  been  redeveloped  based  upon  an  energy 
concept.  These  developments  have  been  guided  by  large-scale  laboratory  tests 
of  disposal  operations  from  a  split-hull  barge  and  a  multiple  bin  disposal  ves¬ 
sel.  Results  from  the  disposal  tests  have  also  been  employed  in  model  valida¬ 
tion  efforts.  Comparisons  of  model  results  with  these  data  are  presented  in  the 
report. 


1  Introduction 


Background 


An  integral  part  of  the  problem  of  managing  a  dredged  material  disposal 
site  is  the  ability  to  determine  the  physical  fate  of  material  immediately  after 
an  individual  disposal  operation  and  ultimately  the  long-term  movement  and/or 
accumulation  of  the  material  deposited  initially  within  the  site.  The  ability  to 
determine  the  short-term  fate  of  dredged  material  disposal  in  open  water  also  is 
an  integral  part  of  assessing  the  water  column  environmental  impact  of 
disposal  operations. 

Field  evaluations  by  Bokuniewicz  et  al.  (1978)  and  laboratory  tests  by 
Johnson,  et  al.,  (1993)  have  shown  that  the  placement  of  dredged  material 
generally  follows  a  three-step  process:  (a)  convective  descent  during  which 
the  material  falls  under  the  influence  of  gravity,  (b)  dynamic  collapse, 
occurring  when  the  descending  cloud  or  jet  either  impacts  the  bottom  or 
arrives  at  a  level  of  neutral  buoyancy,  in  which  case  the  descent  is  retarded 
and  horizontal  spreading  dominates,  and  (c)  passive  transport-dispersion, 
commencing  when  the  material  transport  and  spreading  are  determined  more 
by  ambient  currents  and  turbulence  than  by  the  dynamics  of  the  disposal  oper¬ 
ation.  Mathematical  models  for  predicting  the  short-term  fate  of  material  from 
individual  disposal  operations  that  consider  these  three  phases  have  been 
developed,  e.g.,  Koh  and  Chang  (1973),  Brandsma  and  Divoky  (1976),  and 
Johnson  (1990). 

A  common  deficiency  of  these  numerical  models  is  the  lack  of  data  for 
verification  and  the  inadequacy  of  their  representation  of  the  convective 
descent  and  collapse  phases  in  real  disposal  operations.  For  example,  the 
models  developed  by  Koh  and  Chang  and  subsequently  modified  by  Brandsma 
and  Divoky  and  by  Johnson  treat  the  disposal  from  a  split-hull  barge  as  a 
single  hemispherical  cloud  descending  through  the  water  column. 

As  part  of  the  Dredging  Research  Program  (DRP),  physical  model  disposal 
tests  were  conducted  at  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  (WES).  These  tests  involved  the  disposal  of  various  types  of  material 
from  physical  replicates  of  a  split-hull  barge  and  a  hopper  dredge  in  a  deep 
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basin.  Disposals  were  made  in  water  depths  ranging  from  2.0  to  6.0  ft.^  At  a 
model  scale  of  1:50,  these  tests  simulated  disposal  volumes  of  about  4,000  cu 
yd  from  a  split-hull  barge  and  8,000  cu  yd  from  a  hopper  dredge  in  water 
depths  of  100  to  300  ft. 

Results  from  eight  stationary  and  seven  moving  disposals  from  the  model 
split-hull  barge  and  two  stationary  and  three  moving  disposals  from  the  hopper 
vessel  are  presented  by  Johnson,  et  al.  (1993).  These  results  consist  of  infor¬ 
mation  on  the  short-term  dynamics,  e.g.,  average  descent  and  bottom  surge 
speeds,  suspended  sediment  concentrations  immediately  after  disposal,  and 
bottom  deposition.  Results  from  these  tests  have  been  used  to  guide  the 
numerical  model  developments  discussed  in  this  report  and  to  provide  data  for 
model  verification. 


Purpose 

As  previously  noted,  although  model  development  in  this  area  was  initiated 
in  the  early  1970’s  with  the  work  of  Koh  and  Chang  (1973)  and  was  continued 
with  developments  by  Brandsma  and  Divoky  (1976)  and  Johnson  (1990), 
deficiencies  remained.  Research  in  the  DRP  has  been  directed  at  removing 
many  of  these  model  deficiencies,  e.g.  inadequate  representation  of  disposal 
from  hopper  dredges,  the  inability  to  represent  the  nonhomogenity  of  disposal 
material,  the  inability  to  model  disposals  at  dispersive  sites,  inadequate 
representation  of  the  bottom  collapse  phase,  and  inability  to  model  disposal 
over  bottom  mounds.  With  the  removal  of  these  deficiencies  and  subsequent 
model  verification  employing  the  laboratory  data,  accurate  predictions  of  the 
physical  fate  of  dredged  material  disposed  in  open  water  can  be  obtained. 


Scope 


An  existing  numerical  model  called  DIFID  (Disposal  From  Instantaneous 
Dump)  has  been  extensively  modified  to  yield  a  more  versatile,  accurate  and 
robust  disposal  model  called  STFATE  (Short  Term  FATE).  For  example,  to 
allow  for  disposal  from  hopper  dredges  and/or  a  more  accurate  representation 
of  the  disposal  material,  the  concept  of  multiple  convecting  clouds  that  allows 
for  stripping  of  material  during  descent  has  been  developed.  The  resulting 
model  was  subsequently  verified  using  laboratory  data. 

Part  II  of  the  report  presents  theoretical  details  of  the  STFATE  model.  For 
completeness,  all  theoretical  aspects  of  the  model  are  presented,  although  some 
of  the  discussion  is  a  reproduction  of  discussions  by  Koh  and  Chang  (1973) 
and  Brandsma  and  Divoky  (1976).  Verification  results  from  applying  the 


^  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  (metric)  units  is 
presented  on  page  v. 
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model  to  laboratory  disposal  tests  are  presented  in  Part  III.  Finally,  a  summary 
of  the  model  developments  and  its  verification  along  with  conclusions  from  the 
study  are  given. 
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Regardless  of  the  disposal  method,  the  behavior  of  the  disposal  material  can 
be  separated  into  three  phases:  convective  descent,  during  which  the  disposal 
cloud  or  discharge  jet  falls  under  the  influence  of  gravity;  dynamic  collapse, 
occurring  when  the  descending  cloud  or  jet  either  impacts  the  bottom  or 
arrives  at  a  level  of  neutral  buoyancy  where  descent  is  retarded  and  horizontal 
spreading  dominates;  and  passive  transport-dispersion,  commencing  when  the 
material  transport  and  spreading  are  determined  more  by  ambient  currents  and 
turbulence  than  by  the  dynamics  of  the  disposal  operation.  Figure  1  illustrates 
these  phases. 


Figure  1.  Illustration  of  placement  processes  (from  Moritz  and  Randall 


1992) 

Convective  Descent 

In  STFATE,  multiple  convecting  clouds  that  maintain  a  hemispherical 
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shape  during  convective  descent  are  assumed  to  be  released.  Figure  2  illus¬ 
trates  the  basic  concept. 


By  representing  the  disposal  operation  as  a  sequence  of  convecting  clouds, 
both  split  hull  barge  disposal  as  well  as  disposal  from  a  hopper  dredge  can  be 
modeled.  For  example,  the  material  in  each  hopper  might  be  contained  in  one 
cloud.  This  concept  also  allows  for  a  more  realistic  representation  of  the 
disposal  material  when  consolidation  has  occurred.  More  dense  consolidated 
material  might  be  represented  by  one  cloud  with  the  less  dense  more  fluid-like 
fraction  overlying  the  consolidated  material  represented  by  a  separate  cloud. 

In  addition,  as  illustrated  in  Figure  2,  the  use  of  multiple  convecting  clouds 
also  allows  for  a  more  realistic  representation  of  disposal  from  a  moving 
vessel  where  the  disposal  operation  typically  requires  several  seconds  to 
perhaps  1-2  minutes  for  completion. 

As  each  small  convecting  cloud  descends  through  the  water  column,  mate¬ 
rial  can  be  stripped  and  subsequently  settle  with  its  particle  settling  velocity. 
Movement  of  the  stripped  material  as  small  Gaussian  clouds  is  discussed  later. 

Other  than  the  concept  of  multiple  convecting  clouds  and  the  stripping  of 
material  from  those  clouds,  treatment  of  the  convective  descent  phase  for  each 
of  the  convecting  clouds  is  the  same  as  initially  developed  by  Koh  and  Chang 
(1973).  For  completeness,  the  governing  equations  and  the  following 
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discussion  of  these  equations  are  reproduced  from  Brandsma  and  Divoky 
(1976). 

The  descent  analysis  is  based  upon  the  work  of  Scorer  (1957)  and 
Woodward  (1959)  in  establishing  the  characteristics  of  the  flow  field  in  and 
around  a  buoyant  thermal.  Their  work  treated  clouds  composed  entirely  of 
fluid.  Since  the  solids  concentration  in  discharged  dredged  material  is  usually 
low,  the  cloud  is  expected  to  behave  as  a  dense  liquid,  and  the  buoyant 
thermal  analysis  is  appropriate.  Figure  3  shows  a  schematic  diagram  for  the 
descending  hemispherical  cloud. 


Figure  3.  Convective  descent  (from  Brandsma  and  Divoky  1976) 


A  mean  radius  and  mean  cloud  velocity  are  defined  as  a(t)  and  U  (t) 

where  t  is  time  elapsed  since  release,  p(t)  is  the  mean  density  of  the  element, 
and  pjy)  is  the  ambient  density.  U^(x,y,z,t)  is  the  ambient  current,  which 
is  assumed  to  be  horizontal  and  variable  in  three  dimensions  and  in  time.  The 
characteristics  of  the  cloud  are  assumed  to  remain  similar  at  all  stages  of  the 
convective  descent.  The  various  solid  particles  inside  are  assumed  to  have 
densities  p„  fall  velocities  v^,  and  concentrations  (volume  ratio)  Q.  The 
equations  governing  the  motion  are  those  for  conservation  of  mass, 
momentum,  buoyancy,  solid  particles,  and  vorticity. 

The  time  rate  of  change  of  mass  in  the  cloud  is  the  rate  of  ambient  fluid 
mass  entrainment  minus  the  rate  of  solids  mass  passing  out  of  the  cloud: 
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—  (V  p)  =  Ep  -  YS.  p. 
dt  "  “  7  '  ' 


(1) 


2 

where  the  cloud  volume  is  V  =  _  tt 

3 

The  time  rate  of  change  of  momentum  is  equal  to  the  buoyancy  force 
minus  drag  plus  the  rate  of  ambient  fluid  momentum  entrainment  minus  the 
rate  of  solids  momentum  passing  out  of  the  cloud: 


Several  auxiliary  equations  are  used  for  quantities  in  the  above  equations. 
The  rate  of  entrainment  in  volume  is  the  product  of  the  surface  area  of  the 
hemispherical  front  (2Tra^),  an  entrainment  coefficient  (a),  and  the  magnitude 
of  the  velocity  difference  between  the  cloud  and  the  ambient  fluid: 


E  =  a  \U  -  \ 


(5) 
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The  entrainment  coefficient  will  be  discussed  later  in  conjunction  with  the 
conservation  of  vorticity  equation.  The  volume  rate  of  solids  passing  out  of 
the  cloud  of  the  i‘^  component  is  the  product  of  the  vertically  projected  area  of 
the  cloud,  the  magnitude  of  the  fall  velocity,  the  volume  fraction  of  that 
component  in  the  cloud,  and  a  settling  coefficient: 


Si  =  |Vfi  I 


Several  additional  equations  are  necessary  for  the  following: 


Momentum 


M  =  Cn  1  U 

TYl'  rs 


Buoyancy  force 


F  =  1  -  pj 


Drag  force  in  x-direction 


=  0-5  P 0^^(0-5310^)  ju  -  (u  - 


Drag  force  in  y-direction 


Dy  =  0.5  p^CjjKa^  \U  -  UJ  V 


Drag  force  in  z~direction 


=  0.5  p^C^(0.5jra2)  \U  -  Uj  (w  -  wj 
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Buoyancy 


B  =  ^  ira^  (P«(0)  -  p) 


(12) 


Solid  volume  of  the  i*  component  in  the  cloud 


(13) 


In  the  above  equations,  a  is  an  entrainment  coefficient;  )8,  is  a  settling 
coefficient;  j  is  the  unit  vector  in  the  vertical  direction;  Q,  is  a  drag 
coefficient;  C„  is  an  apparent  mass  coefficient;  and  pjo)  is  the  density  at  the 
free  surface.  The  drag  coefficient  Q,  is  a  function  of  Reynolds  number  and 

therefore  depends  on  a  and  \IJ  -  U^\.  The  value  of  Q,  is  suggested  to  be 
0.5.  The  apparent  mass  coefficient,  C„,  is  suggested  to  be  between  1.0  and 
the  value  for  a  sphere,  i.e.,  1.5. 

The  last  governing  equation  is  that  for  vorticity.  The  total  vorticity  is  the 
cloud’s  identity-preserving  mechanism,  and  it  is  also  important  in  determining 
the  amount  of  entrainment.  When  a  cloud  of  material  is  ejected  into  the 
ambient  fluid,  some  initial  vorticity  is  generated  when  passing  through  the 
entrance  boundary.  Total  vorticity  is  generated  only  by  shear  forces  at  the 
fluid  boundaries.  Once  a  cloud  is  in  the  ambient  fluid,  there  are  two  pos¬ 
sibilities  if  the  bottom  is  not  encountered.  In  a  uniform  density  ambient  fluid, 
total  vorticity  is  conserved,  although  cloud  growth  acts  to  diffuse  the  vortex 
strength.  In  a  stratified  fluid,  the  density  gradient  acts  to  decay  the  total 
vorticity  according  to: 


dK 

dt 


-  A  e 


(14) 


where  K  is  vorticity;  y4  is  a  dissipation  parameter;  and  e  is  the  density  gradient 
as  defined  below: 


a  =  £^ 

P.(0) 


(15) 
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e 


dp  a 
~df 


(16) 


C  is  a  vorticity  dissipation  coefficient  which  is  equal  to  3  according  to  Turner 
(1960).  Vorticity  is  of  concern  only  because  it  influences  entrainment. 

The  entrainment  coefficient  a  is  dependent  upon  the  properties  of  the 
cloud,  the  properties  of  the  ambient  fluid,  and  the  turbulence  structure  inside 
and  outside  the  convecting  cloud.  In  formulating  an  expression  for  an 
entrainment  coefficient,  it  is  necessary  to  account  for  the  structure  of  the  cloud 
as  it  changes  from  a  vortex  ring  to  a  turbulent  thermal.  Scorer  (1957)  and 
Richards  (1961)  experimentally  determined  the  entrainment  coefficient  for 
turbulent  thermals,  a^,  to  be  approximately  0.25.  In  studying  the  motion  of  a 
vortex  ring.  Turner  (1960)  found  an  entrainment  coefficient: 


B 

a  =  - 

23t  gC^K^ 


(17) 


by  assuming  similarity  where  Cj  was  found  to  be  0.16.  B  is  the  buoyancy  and 
K  is  the  total  vorticity.  As  the  cloud  descends  and  its  vorticity  approaches 
zero,  Turner’s  assumption  of  similarity  cannot  hold.  Since  a  is  expected  to 
approach  found  in  turbulent  thermals,  Koh  and  Chang  (1973)  thought  it 
reasonable  to  postulate  that  the  dependence  of  a  on  5  and  K  might  be  of  the 
form: 


a 


a. 


N 


tanh 


B 


2KgC^K'^aQ 


(18) 


Their  only  justification  for  Equation  18  was  that  it  tends  to  the  correct  limits; 
to  Turner’s  relation  when  K  is  large  and  to  when  K  is  small. 

Laboratory  studies  by  Bowers  and  Goldenblatt  (1978)  resulted  in  the  fol¬ 
lowing  empirically  derived  expression  for  the  entrainment  coefficient. 


a  =  0.285  +  0.00493  {MIL  -  2.9) 


(19) 


where  MLL  =  Multiple  of  Liquid  Limit.  This  expression  results  in  a  maxi¬ 
mum  value  of  about  0.30.  However,  Kraus  (1991)  determined  entrainment 
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coefficients  of  about  0.60  from  field  data  collected  at  Mobile  Bay.  As  noted 
later,  the  larger  value  is  required  to  achieve  volume  increases  of  50  to  70  times 
the  disposal  volume  as  observed  by  Bokuniewicz,  et  al.  (1978)  and  Johnson, 
et  al.  (1993). 

Koh  and  Chang  used  dimensional  analysis  to  show  that  the  dimensionless 
mass  rate  of  settling  is  a  function  of  the  ratio  of  the  descending  velocity  of  the 
cloud,  v;  the  fall  velocity  of  the  solid  particles,  v^;  the  concentration  of  each 
group  of  particles,  C^-;  and  total  concentration,  C: 

_J_  =  JiC,,(l-P,)  (20) 


where  is  defined  to  be  a  settling  coefficient  which  depends  upon  v/v^j ,  C^-, 
and  C.  P,-  is  expected  to  be  between  0  and  1,  representing  the  two  cases  of 
settling  freely  or  no  settling. 

Therefore, 


Pf  = 


b  if 


if 


<  1 


&  1 


(21) 


Given  a  set  of  initial  conditions,  e.g.,  the  total  volume  of  disposal  material, 
sediment  concentrations,  etc..  Equations  1-4  along  with  Equation  14  can  be 
solved  using  numerical  techniques.  In  STFATE,  the  forward  Euler  scheme  is 
used.  Earlier  work  by  Koh  and  Chang  employed  a  fourth  order  Runga-Kutta 
scheme.  However,  testing  over  a  range  of  disposal  volumes  and  water  depths 
revealed  little  difference  in  the  computed  results  and  the  first  order  Euler 
scheme  appeared  to  be  more  robust.  Given  the  idealized  representation  of  the 
disposal  operation  and  the  uncertainty  of  initial  conditions,  the  need  for  higher- 
order  numerical  schemes  is  questionable. 

Bokuniewicz,  et  al.  (1978)  found  that  the  insertion  speed  was  an  important 
variable.  Therefore,  special  attention  has  been  directed  at  computing  the  cloud 
insertion  speed  in  STFATE. 

Consider  Figure  4,  with  being  the  inside  area  of  the  disposal  vessel, 
being  the  area  of  the  opening,  h  the  height  of  the  disposal  material,  d  the  ves¬ 
sel  draft,  the  atmospheric  pressure,  and  and  Vj  velocities  at  the  surface 
and  at  the  opening  respectively. 
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Figure  4.  Computation  of  insertion  speed 


The  energy  equation  can  be  written  as 


(22) 


where,  at  the  moment  the  doors  are  opened,  the  pressure  at  the  opening,  is 


^1=^0  + 


(23) 


Therefore, 
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(24) 


Tr2  X  r2 

h  +  _  =  _  d  +  _  +  e 

2g  P  *  2g 


where  e  is  the  energy  lost  due  to  friction  as  the  material  passes  through  the 
opening. 

From  continuity, 


Thus  Equation  24  can  be  solved  for  Vj  to  yield 


2g  h-^d-z 

P*  j 


1  -  _L 
A„ 


Before  the  doors  are  opened,  the  mass  of  the  boat  and  the  disposal  material 
must  equal  the  mass  of  the  water  displaced,  i.e., 

dp  ^  =  hp^g  +  mass  of  boat 


P  w  ^ 

h  =  -  unloaded  draft 

Ps^o 


where  A  =  area  of  the  vessel  at  the  water  line.  Therefore,  as  convecting  clouds 
are  released,  changes  in  h  and  d  are  related  by 
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where  is  the  average  or  bulk  density  of  the  entire  load,  Volj  is  the  total 
volume  being  disposed,  and  and  Vol’  are  the  average  density  and  volume  of 
the  individual  layers  in  the  disposal  vessel.  Each  layer  ultimately  becomes  a 
convecting  cloud. 

Dynamic  Collapse  in  Water  Column 

As  the  disposal  cloud  goes  through  the  convective  descent  phase,  it  gains 
mass  and  momentum  by  entrainment.  The  horizontal  velocity  of  the  cloud 
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tends  to  approach  that  of  the  ambient  fluid.  Coincidently,  the  waste  material 
concentration  is  greatly  reduced  and  the  vorticity  becomes  insignificant 
because  of  dissipation  by  ambient  stratification  and  turbulence.  If  the  cloud 
reaches  the  depth  of  neutral  buoyancy,  its  momentum  will  tend  to  make  it 
overshoot  beyond  the  neutrally  buoyant  point  while  the  buoyancy  force  will 
tend  to  bring  it  back  to  the  neutrally  buoyant  position.  The  combined  action 
of  these  forces  will  make  the  cloud  undergo  a  decaying  vertical  oscillation.  As 
the  vertical  motion  of  the  cloud  is  being  suppressed,  the  cloud  tends  to 
collapse  vertically  and  spread  out  horizontally,  seeking  hydrostatic  equilibrium 
within  the  stratified  ambient  fluid.  As  the  cloud  collapses,  its  cross  section 
becomes  elongated  in  the  horizontal,  and  another  dimension  is  needed  to 
describe  the  cloud  shape.  If  the  cross  section  of  the  cloud  is  assumed  to  be 
elliptical,  its  shape  may  be  characterized  by  its  semi-major  and  semi-minor 
axes,  b  and  a,  respectively  (Figure  5). 

It  is  assumed  that  the  cloud  always  retains  the  shape  of  an  oblate  spheroid. 
If  coordinate  axes  are  chosen  to  originate  from  the  cloud  centroid,  the  cross- 
sectional  outline  of  the  cloud  may  be  represented  by 


y 


a 


a 


2 


(33) 


where  a  and  h  vary  with  time.  The  cloud  is  assumed  to  remain  symmetrical, 
which  in  practice  will  only  be  true  if  there  is  no  relative  velocity  between  the 
cloud  and  the  ambient  fluid  and  therefore  no  velocity  shear.  Following  con¬ 
vective  descent,  the  velocity  difference  between  the  cloud  and  the  ambient 
fluid  is  expected  to  be  small,  and  its  influence  on  shape  can  be  neglected. 

With  the  exception  of  vorticity,  the  conservation  equations  used  for  con¬ 
vective  descent  still  hold.  Any  differences  are  due  to  the  additional  dimension 
used  to  describe  the  cloud.  The  conservation  equations  are: 

Conservation  of  mass 


l(V,p)=£p„-K,.p,. 


(34) 


where  the  cloud  Volume  V^,  = 


—  n  ab^. 
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Figure  5.  Aspects  of  cloud  collapse  (from  Brandsma  and  Divoky  (1976)) 


Conservation  of  momentum 


^  D  +  C/7-  ES.  p  .  C/ 

dt  i 


(35) 


The  major  auxiliary  equations  are  those  for  entrainment  of  ambient  fluid, 
for  the  settling  of  solid  particles,  and  for  the  collapse  of  the  cloud. 

The  two  contributions  to  entrainment  come  from  movement  of  the  entire 
cloud  through  the  ambient  fluid  and  from  the  additional  velocity  shear  at  the 
cloud  boundary  due  to  the  cloud  collapse.  Each  contribution  is  the  product  of 
the  surface  area  of  the  cloud,  an  entrainment  coefficient,  and  a  velocity.  The 
expression  for  the  rate  of  entrainment  of  ambient  fluid  is: 


where  R  =  yb^-a^  ;  a  is  the  entrainment  coefficient  for  cloud  motion;  and 

I  U  -  U  g  I  is  the  magnitude  of  the  velocity  difference  between  the  cloud  and 
the  ambient  fluid,  is  introduced  into  the  formulation  to  account  for  the 
entrainment  due  to  collapse,  and  db/dt  is  the  velocity  of  the  tip  of  the  collaps¬ 
ing  cloud.  It  is  assumed  that  entrainment  due  to  the  cloud  motion  should  die 
out  as  the  cloud  settles  into  the  neutral  buoyancy  position.  This  may  be  done 
by  letting  the  entrainment  coefficient  be: 


Chapter  2  Model  Theory 


17 


a  = 


(39) 


where  olq  is  the  entrainment  coefficient  for  a  turbulent  thermal.  This  relation 
is  employed  merely  as  an  effective  way  to  turn  off  entrainment  due  to  cloud 
motion. 

As  in  the  section  on  convective  descent,  the  volume  rate  of  solids  settling 
out  of  the  cloud  is 


5,  =  |v,|  C,  (1  -  ^,) 


(40) 


where  |(3,  is  the  settling  coefficient  defined  by  Equation  21. 

The  mechanism  that  drives  the  collapse  of  the  cloud  within  the  water 
column  is  the  density  difference  between  the  inside  and  outside  of  the  cloud. 

It  is  assumed  that,  because  of  the  turbulent  mixing,  the  density  gradient  inside 
the  cloud  is  less  than  that  outside.  Assume  that  the  cloud  is  resting  at  the 
level  of  neutral  buoyancy  and  that  the  ambient  density  at  this  level  is  Pq. 
Assume  further  that  the  density  gradients  inside  and  outside  of  the  cloud  are 
of  constant  magnitude.  Let  e  be  the  normalized  ambient  density  gradient: 


Po  ^ 


The  density  gradient  inside  the  cloud  is  assumed  to  be  less  than  that  in  the 
ambient  fluid  by  a  factor  yag  la,  where  7  is  a  coefficient,  Cg  is  half  the  final 
radius  of  the  convective  descent  cloud  (Figure  5),  and  a  is  the  semi-minor  axis 
of  the  collapsing  cloud.  Following  these  assumptions,  the  ambient  density  in 
the  region  of  the  cloud  is 

P^- Pg{l  -  f:y')  =  Pg{\-e{a-y))  (42) 


and  the  density  inside  the  cloud  is 
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P  =  Po 


(43) 


\ 

1  -  Y  _  ey' 

=  Po 

1  -Y  —  e(a-y) 

"  J 

where  y  and  y’  are  as  defined  in  Figure  6.  These  density  profiles  are 
illustrated  in  Figure  5b. 

The  centroid  of  the  collapsing  slice  of  the  cloud  shown  in  Figure  5c  moves 
with  respect  to  the  centroid  of  the  cloud.  Consider  the  slice  (of  angular 
dimension  dQ)  as  a  free  body  and  integrate  the  pressures  over  the  surfaces  of 
the  slice  to  obtain  the  radial  force  (acting  at  the  slice  centroid)  driving  the 
collapse  of  the  slice.  The  pressures  are  assumed  to  be  hydrostatic.  The  pres¬ 
sure  on  the  rounded  external  surface  of  the  slice  is  simply  the  pressure  in  the 
ambient  fluid,  pjiy),  and  it  is  integrated  over  the  projected  area  as  shown  in 
Figure  6. 


^ext  =  Jp*  Pa^)  '■(y)  ^0 


(44) 


Figure  6.  Slice  of  collapsing  cloud,  showing  integration  elements  for 

determining  the  driving  force  of  collapse  (from  Brandsma  and 
Divoky  (1976)) 

The  pressure  inside  the  cloud,  Pc(y,r),  is  a  function  of  the  vertical  and  radial 
position  inside  the  cloud.  This  pressure  is  integrated  over  both  side  faces  of 
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the  slice,  and  dQ  is  assumed  small  to  give  an  expression  for  the  radially 
outward  directed  force: 


(45) 


where  R  = 


Once  the  integrals  (44)  and  (45)  are  evaluated  and  the  difference  of 
is  taken,  the  radial  force  driving  collapse  is  obtained: 


^0  3 

JtPo(l-Y— )  b 

Fjy  =  - f -  dQ 

"  16 


(46) 


1  3p 

Since  e  =  _  _  ,  the  expression  may  be  written  as 

Po 


(i^)  gaH  d6 

16  a  dy 


(47) 


In  formulating  the  inertia  of  the  slice  and  force-resisting  collapse,  it  is 
assumed  that  the  horizontal  velocities  of  the  elements  inside  the  slice  are 
related  to  the  radial  distance,  r,  from  the  centroid  of  the  cloud  and  that  the 
velocity  of  horizontal  deformation  is  characterized  by  the  velocity  of  the 
centroid  of  the  slice  segment.  The  cloud  tip  velocity  due  to  collapse  (which 
will  be  called  v^)  is  linearly  related  to  the  slice  centroid  horizontal  velocity  by 


16  dc 
3n  dt 


(48) 


The  forces  resisting  the  collapse  of  the  slice  are  form  drag,  Dj^,  and  skin  fric¬ 
tion  drag  F^ 
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(49) 


Dd  *  Pa  Kll  Vi  dB 


(50) 


where  Q  is  the  drag  coefficient  for  a  wedge  and  Cj^i^  is  the  friction  coeffi¬ 
cient  for  a  flat  plate,  which  is  a  function  of  the  kinematic  viscosity  (Hoemer, 
1965). 

The  horizontal  inertia  of  the  cloud  slice  is  the  time  rate  of  change  of  the 
product  of  its  mass  and  the  velocity  of  the  slice  centroid; 


I  = 


nab^ 

“TT" 


dQ 


(51) 


The  dynamic  equation  is  formulated  as  the  summation  of  the  external  forces 
acting  on  the  slice  which  is  equated  to  the  inertia. 

I  ^Fn-DD-Ff  (52) 


The  integral  of  Equation  52  over  one  time  step  is  used  to  determine  the  value 
of  v^  at  the  end  of  the  time  step  with  d©  canceling  out  of  all  terms. 

The  effects  of  entrainment  on  mass,  momentum,  and  buoyancy  have  been 
considered  previously.  Entrainment  also  influences  the  shape  of  the  cloud 
cross  section,  which  is  assumed  to  remain  symmetric.  The  rate  of  volume 
increase  of  the  cloud  due  to  entrainment  is  assumed  to  be  accounted  for  by  an 
increment  in  the  tip  velocity  with  the  vertical  dimension  of  the  cloud  held 
constant.  This  increment  will  be  called  ^2-  'Th®  velocity  of  the  slice  tip  of  the 
collapsing  and  entraining  cloud  is: 


db 

dt 


+  V2 


(53) 


Equation  (34),as  it  stands,  accounts  only  for  the  gross  amount  of  mass 
entrained.  It  may  be  rewritten  under  the  assumptions  of  constant  cloud  density 
and  constant  vertical  dimension  during  each  integration  step  to  describe  how 
entrainment  adds  to  the  growth  of  b: 
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(54) 


£p,-a,Pi 


"2  = 


p  ab 
3 


Several  additional  equations  are  used  for: 
Momentum 


Cj^p^nab'^tf 


(55) 


Buoyancy  force 

F  =  _p  j 


(56) 


Drag  force  in  the  x-direction 


1 

2 


p^Co^Jtab  \U  -U^\  (m-mJ 


(57) 


Drag  force  in  the  y-direction 


\U  -  U^\  V 


(58) 


Drag  force  in  the  z-direction 


p^CD^ab'^  ICT-^I  (»v-H.J 


(59) 
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Buoyancy 


5  =  ljtaft2(p^(0)-p)  (60) 

Solid  volume  of  the  component 

Pi  =  (^1) 


In  the  above  equations,  Cj^  is  an  apparent  mass  coefficient,  Cjy^  is  the  drag 
coefficient  for  a  spheroid^  wedge,  and  is  the  drag  coefficient  for  a 
circular  plate. 

As  in  the  convective  descent  phase.  Equations  34-37  are  solved  by  the 
forward  Euler  scheme  to  yield  the  time  history  of  the  disposal  cloud  while 
collapsing  in  the  water  column.  Initial  conditions  such  as  the  bulk  density  and 
axes  of  the  oblate  spheroid  are  provided  from  the  last  time  step  of  the 
convective  descent  phase. 


Dynamic  Collapse  on  Bottom 

As  previously  noted,  the  convective  descent  and  water  column  collapse 
phases  are  essentially  the  same  as  developed  by  Koh  and  Chang  (1973)  and 
virtually  all  of  the  discussion  presented,  with  the  exception  of  the  discussion 
on  insertion  speed,  has  been  taken  from  Brandsma  and  Divoky  (1976).  How¬ 
ever,  the  manner  in  which  collapse  of  the  disposal  cloud  is  treated  when  the 
bottom  is  encountered  during  convective  descent  is  completely  different. 

Bottom  collapse  in  STFATE  is  computed  from  a  conservation  of  energy 
concept  When  the  cloud  strikes  the  bottom,  it  possesses  a  certain  amount  of 
potential  energy  which  can  be  computed  since  the  mass  of  the  cloud  and  the 
location  of  its  centroid  are  known.  In  addition,  the  kinetic  energy  of  the 
impacting  cloud  can  be  computed  since  its  velocity  is  known.  Thus,  the  total 
energy  of  the  cloud  at  the  moment  of  impact  is  known.  This  energy  is  then 
available  to  drive  the  resulting  bottom  collapse  or  surge. 

A  basic  assumption  is  that  the  bottom  collapsing  cloud  is  one-half  of  an 
ellipsoid.  If  the  bottom  is  flat,  the  cloud  becomes  one-half  of  an  oblate 
spheroid.  From  Figure  7,  the  general  equation  of  an  ellipsoid  is 
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Figure  7.  Bottom  collapse 


X 


2 

+ 


(62) 


The  kinetic  energy  of  the  collapsing  cloud  can  be  written  as 


KE  = 


i  +  w^)  dV 


(63) 


where  p  is  the  cloud  density,  (m,v,m')  represent  the  speed  of  any  part  of  the 
cloud  along  the  (jc,y,z)  axes,  respectively,  and  V  is  the  volume  of  the  cloud. 
Expressions  for  (m,v,w)  can  be  deduced  by  integrating  Equation  62  with 
respect  to  time,  yielding 


(64) 
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Since,  by  definition,  u  =  x  ,  v  =  y  ,  and  w  =  z  ,  for  Equation  64  to  be 
satisfied  the  following  holds: 


(66) 


Substituting  Equations  65-67  into  Equation  63  yields  the  following  expression 
for  the  kinetic  energy: 


Now  consider  the  evaluation  of  the  general  expression 

x'-dV-  j‘x-Adx 
V 

where  A  is  the  cross-section  area  in  the  yz  plane.  Considering  cross-sections 
in  the  xz  and  xy  planes,  the  following  equations  for  an  ellipse  can  be  written: 
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In  addition,  the  area  can  be  written  as 


A  = 


(72) 


Solving  Equation  70  for  z  and  Equation  71  for  j  and  substituting  into 
Equation  72,  yields 


A 


(73) 


Substituting  Equation  73  into  Equation  69  and  evaluating  the  resulting 
integrals  yields 


x"  dV  =  1  acb"*' 

4 


1 


1 

n+3 


(74) 


Using  Equation  74  to  evaluate  the  integrals  in  Equation  68  and  now  letting  V 

=  —  Trabc  be  the  volume  of  a  half  ellipsoid,  the  kinetic  energy  for  a 

3 

quadrant  of  the  half  ellipsoid  becomes 


1  -2  -2  -2, 
KE  =  -L  pV  (b  +  a  +  c) 


(75) 


An  expression  for  the  potential  energy  can  be  derived  as  follows.  Consid¬ 
ering  Figure  8,  the  total  potential  energy  of  the  surrounding  ambient  fluid  and 
the  collapsing  cloud  is 
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Figure  8.  Potential  energy  of  bottom  cloud 


where  Vj-  =  Afl  =  {A^  -  A)H  +  AH 
Equation  (76)  can  now  be  written  as 

=  1 1 1  pgy  dV^  +  pgy  d\ 

{aI-a)h  ah 

The  second  integral  in  Equation  77  can  be  written  as 
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■  [  pgydv,  =  f  f 

a  H  ~ 

pg  {ydy  +  p^g  [  ydy 

J  J  J 

A  r 

V  a 

^ 

I  {p-p^8  I  ydy 


I  Pas\ydy 


dA 


dA 


dA 


(78) 


or 


AH 


PaSy 


dV.  + 


A.S|J| 

Aa 


dV 


(79) 


If  the  free  surface  isn’t  moving,  the  first  integral  above  is  a  constant.  Equa¬ 
tion  74  can  be  used  to  evaluate  the  second  integral,  yielding  the  following 
expression  for  the  potential  energy  of  a  quadrant  of  the  half  ellipsoid. 


PE  =  —  ApgVa  +  constant 


(80) 


Since  only  changes  in  potential  energy  influence  the  dynamics  of  the 
collapsing  cloud,  the  constant  does  not  need  to  be  evaluated. 

Given  that  the  change  in  the  total  energy  of  the  cloud,  which  is  the  sum  of 
the  kinetic  and  potential  energies,  is  equal  to  work  done  by  the  cloud,  the 
following  expression  holds: 


AKE  +  APE  =  WORK 


(81) 


What  is  desired  are  expressions  for  the  time  rate  of  change  of  a,  b  and  c. 
Assuming  that 
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(82) 


and  substituting  into  Equation  75,  yields 


XE  =  _pV 


{bf  +  {af  + 


{bf 


(83) 


Since  the  time  rate  of  change  of  the  volume  of  the  half  ellipsoid  is  the  rate  of 
entrainment  of  ambient  fluid,  Q^,  the  following  expression  can  be  derived: 


a  =  a 


Qe 

V 


b 

b 


(84) 


Substituting  Equation  84  into  Equation  83  yields  the  following  equation  from 
which  the  change  of  b  each  time  step  can  be  determined. 


2  -2  2a2  a'^b^^ 


,  b^  a  ^ 
b^ 


'  Ab' 

At 

22. 

r  1  b] 

'  Ab' 

V 

b 

At 

2^2 


\QKE  «  2. 


(85) 


pV 


V' 


Once  Afe  is  determined,  the  change  in  c  is  determined  from  Equation  82  and 
the  change  in  a  from  Equation  84. 

The  computations  proceed  in  the  following  manner.  At  the  end  of  the 
convective  descent  phase,  the  kinetic  and  potential  energies  of  the  cloud  are 


KE  =  PttR^  \U\^ 


(86) 
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=  I  ApgR^a 


(87) 


where  R  is  the  radius  of  the  hemispherical  cloud  and  a=b=c=R  at  the  initia¬ 
tion  of  collapse.  The  exception  is  when  disposal  takes  place  in  shallow  water. 
In  that  case,  a  is  set  to  be  the  distance  from  the  bottom  of  the  disposal  vessel 
to  the  sea  bottom  and  &=c  is  determined  so  as  to  conserve  the  volume  of 
disposal  material. 

The  potential  energy  of  the  cloud  changes  as  the  cloud  density  decreases 
due  to  entrainment  of  ambient  fluid  and  the  loss  of  sediment,  as  well  as,  due 
to  the  height  of  the  centroid  of  the  cloud  decreasing  as  the  cloud  spreads  over 
the  bottom.  The  change  in  kinetic  energy  is  computed  from  Equation  81  once 
the  work  done  by  the  collapsing  cloud  is  quantified. 

Work  must  be  done  to  overcome  bottom  friction,  drag,  the  production  of 
internal  mrbulence  and  setting  the  ambient  fluid  in  motion.  The  contribution 
to  the  work  term  in  Equation  81  due  to  bottom  friction,  Wp  is 


dt 


dc 

dt 


dt 


dc 

dt 


At 


The  contribution  due  to  drag,  W^,  is 


=  \pa  Cjras 


- 

[dbV  ^ 

I  J 

'  dcV 

dt 

\ 

[dbV  ^ 

1  dt\ 

'  dcV 
dc 

(88) 


(89) 


The  work  required  to  set  the  ambient  fluid  in  motion  is  assumed  to  be  some 
percent  of  the  existing  kinetic  energy,  with  the  value  being  input  as  a  model 
coefficient.  The  contribution  due  to  the  generation  of  internal  turbulence,  W„ 
is 
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(90) 


W,  =  —  Tt  abc  p  u*^  —  At 
‘3  ^  dz 

where  the  friction  velocity,  u*,  is  computed  from 


where  is  the  height  of  the  bottom  roughness.  The  velocity  gradient  in 
Equation  90  is  computed  from 


dz  a 


Once  the  change  in  kinetic  energy  is  computed,  the  new  kinetic  energy  is 
known  and  the  changes  in  c,  a,  and  b  can  be  determined  from  Equations  82, 
84  and  85,  respectively. 

Time  Steps 

The  computational  time  step  in  both  the  convective  descent  and  collapse 
phases  is  determined  internally  within  STFATE.  For  the  case  of  a  uniform 
ambient  density,  the  first  attempt  at  setting  the  time  step  in  the  convective 
descent  phase  is 


whereas,  if  the  ambient  density  is  non-uniform 
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0.01  Jt  w 


DT 


1 

/ 

\ 

^9  a  R 

P  J 

v 

J  Ap  //  J 

(94) 


If  the  total  number  of  time  steps,  ISTEP,  in  the  convective  descent  phase  is 
less  than  100  or  more  than  200,  the  computations  are  reinitialized  with  the  new 
time  step  determined  from 


DT. 


{DT)  (ISTEP) 
150 


(95) 


Up  to  ten  trials  are  allowed.  Extensive  testing  of  STFATE  over  a  wide  range 
of  conditions  have  shown  that  the  computations  usually  converge  within  ten 
trials.  If  not,  results  from  the  10***  trial  are  used. 

If  collapse  occurs  in  the  water  column,  the  initial  time  step  is  computed 
from 


DT 


0.001 


840 

P  dy 

0.42857' 

R 

0.1 


N 


g_  ^ 
P  3y 


(96) 


whereas,  if  collapse  occurs  on  the  bottom,  it  is  set  to  be  that  used  in  the  con¬ 
vective  descent  computations.  If  the  total  number  of  collapse  time  steps, 
ISTEPC,  is  less  than  100  or  greater  than  400,  the  computations  are  reinitialized 
and  the  new  time  step  is  computed  from 


DT. 


(DT)  (ISTEPC) 
250 


(97) 


Equations  93,  94  and  96  were  employed  by  Koh  and  Chang  (1973). 
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Bottom  Slope 


The  influence  of  a  bottom  slope  is  to  increase  the  spread  of  the  bottom 
surge  in  a  downslope  direction  and  to  decrease  the  spread  in  an  upslope  direc¬ 
tion.  Within  the  idealized  framework  of  a  collapsing  half  ellipsoid,  this  effect 
is  modeled  in  the  following  manner  in  STFATE. 

Consider  Figure  9.  Slopes  sxl,  sx2,  szl  sz2  are  computed  at  the  four  cen¬ 
troids  shown  by  subtracting  the  depth  at  the  cloud  centroid  from  the  depth  at 
the  individual  centroids.  Movement  of  the  four  centroids  shown  due  to  a  slope 
can  then  be  determined  from 


Figure  9.  Effect  of  bottom  slope 


— 1  =  g  sin  (sxl)  (98) 

dt  p 

— i  =  g  sin  (sx2)  (99) 

dt  p 

— 1  =  g  sin  (szl)  (100) 

dt  p 
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Once  Up  Up  Wj,  and  W2  are  determined,  the  change  in  bp  bp  Cp  and  C2  can  be 
determined  from 


> 

11 

8 

(102) 

At 

'  3 

1 

Ab2 

8 

(103) 

“aT 

'  T 

2 

Ac  2 

8 

(104) 

"aT  ' 

'  T 

1 

Ac2  _ 

8 

(105) 

“aT  ' 

3 

z 

The  contribution  to  the  bottom  spreading  as  a  result  of  bottom  slope  then 
becomes 


Ab  =  Abj  +  ^2 

(106) 

Ac  =  ACj  +  Ac2 

(107) 

Since  the  time  rate  of  movement  of  each  of  the  four  centroids  is  equal  to  the 
speeds  determined  in  Equations  98-101,  the  new  locations  of  the  four  centroids 
can  be  determined.  The  new  location  of  the  total  cloud  centroid  is  then  deter¬ 
mined  as  an  average  of  the  locations  of  the  individual  centroids. 


Disposal  at  Dispersive  Sites 

In  the  earlier  disposal  model  work,  there  was  no  mechanism  for  keeping 
material  from  depositing,  and  once  material  was  deposited  on  the  bottom  it 
remained  there.  Thus,  at  dispersive  sites,  model  results  were  questionable.  In 
STFATE,  a  critical  shear  stress  for  deposition  is  specified  for  each  sediment 
fraction.  A  computed  bottom  shear  stress  is  then  compared  to  the  critical 
value  to  determine  if  suspended  material  can  be  deposited. 
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If  the  bottom  shear  stress  is  affected  only  by  currents  and  no  waves  are 
present,  it  is  computed  from  the  following  expression  based  upon  the 
assumption  of  a  log  profile. 


=  Pa 


0.4  ^ 


2  2 


in 


I  \ 
z 


2 

=  Pa»* 


(108) 


where  (u^  w^)  are  the  components  of  the  ambient  velocity,  z  is  the  height  at 
which  (u^  wjl  are  determined,  e.g.  1  m,  and  is  the  bottom  roughness,  e.g.,  1 
mm.  M*  is  the  friction  velocity  and  becomes  if  waves  are  present. 

If  waves  are  present  at  the  disposal  site,  the  following  approach  taken  from 
Madsen  and  Wikramanayake  (1991)  is  used.  With  this  approach,  the  ambient 
current  angle,  wave  amplitude,  wave  period,  and  an  initial  value  for  the  wave- 
current  angle,  wc^,  must  be  input.  The  wave-current  friction  factor,  is 
computed  first  from 


^  +  logjo 

1  ) 

=  logic 

4  \ 

30z 

+  0.1 


(109) 


if  A^/(30  zj  >  1000,  where  is  the  wave  amplitude.  However,  if  A^,/(30zJ 
<  1000,  the  following  expression  is  employed. 


/wc  =  exp 


5.2 


30z. 


-0.19 


-  6.1 


0.24 


\-1.2 


302. 


(110) 


The  next  step  involves  an  iteration  for  the  current  friction  velocity  and  the 
wave-current  angle.  The  initial  guess  for  the  current  friction  velocity  is  m*^. 
As  noted,  the  initial  guess  for  the  current-wave  angle  is  input.  The  followmg 
iteration  involving  Equations  113  and  116  is  then  performed  and  continues 
until  the  values  of  the  current  friction  velocity  and  wave-current  angle,  wc, 
converge,  i.e.. 


±  >  0.99 


’»-i 


(111) 
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(112) 


wc 


wc 


i-1 


>  0.99 


where  i  is  the  iteration  number.  The  procedure  is  as  follows.  First  compute  II 
from 


/  N 

{  2  u,  \ 

*/-l 

4 

in 

0.4 
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\  /  1 
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(113) 


where 
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wc 

IJQ 


(114) 


C  = 
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o 


30 
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*  w 


(115) 


and  is  the  input  wave  period.  Then  compute  12  from 


12 


'i-l 


0.4 


In 


0.8 

T 


0.53  M, 

(1-0  -C)«. 

-  0.78  +  a 


\ 


/ 


(116) 


The  new  estimates  for  the  current  friction  velocity,  u,.,  and  the  wave-current 
angle,  wCj,  at  the  i-th  iteration  are  then  computed  from 
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(117) 


u 


^-1 


/ 

\ 


where 


ir 


(118) 


and 

wc^  =  +  (wc^  “  a) 


(119) 


where 


a  =  TAN-'^ 


II 

im 


(120) 


With  the  convergence  of  m*.  and  wc^-,  the  total  bottom  shear  stress  is  then 
computed  from 


1/2 


cos^wc^  +  p  ^  sin^H'Cjj 


(121) 


The  shear  stress  computed  from  either  Equation  108  or  Equation  121  is  then 
compared  to  the  critical  shear  stress  for  deposition  to  determine  if  a  particular 
sediment  fraction  deposits. 


Transport-Diffusion  Phase 

At  most  disposal  sites,  the  convective  descent  and  dynamic  collapse  phases 
only  last  on  the  order  of  a  few  minutes.  When  the  rate  of  spreading  of  the 
collapsing  cloud  becomes  less  than  an  estimated  rate  of  spreading  due  to 
turbulent  diffusion,  the  collapse  phase  is  terminated  and  the  "longer"  term 
transport-diffusion  phase  is  initiated.  In  this  phase,  material  in  suspension  is 
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transported  and  diffused  by  the  ambient  current  while  undergoing  settling.  Of 
course  any  non-sediment  constituent  being  modeled  is  also  transported  and 
diffused. 

A  basic  assumption  in  the  transport-diffusion  phase  is  that  material  (or 
constituent)  concentrations  can  be  determined  from  a  superposition  of  small 
clouds  characterized  by  a  normal  or  Gaussian  distribution,  i.e., 


C 


m 


exp 


-1/2 


2 


where 

m  =  volume  of  solids  in  the  cloud, 

=  standard  deviations,  ft 
x,y,z  =  spatial  coordinates,  ft 
Xo.yo.Zo  =  coordinates  of  cloud  centroid,  ft 


(122) 


These  clouds  are  formed  as  material  is  stripped  away  during  the  descent  of  the 
convecting  clouds  as  well  as  during  the  collapse  phase.  For  collapse  in  the 
water  column,  small  clouds  are  formed  as  material  settles  from  the  collapsing 
cloud.  However,  during  collapse  on  the  bottom,  laboratory  experiments  by 
Johnson,  et  al.  (1993)  as  well  as  field  data  collected  at  Mobile,  AL  by  Kraus 
(1991)  imply  that  fine  material  is  also  lost  to  the  water  column  at  the  top  of 
the  collapsing  cloud. 


At  the  end  of  each  time-step,  each  cloud  is  advected  horizontally  by  the 
input  velocity  field.  The  new  position  of  the  cloud  centroid  is  determined  by 


X  =  X  +  u  •  At  (123) 

z„  =  z„  +  w  •  At  (124) 


where 


u,w  =  input  ambient  velocities,  fps 
At  =  long-term  time-step,  sec 

In  addition  to  the  advection  or  transport  of  the  cloud,  the  cloud  grows  both 
horizontally  and  vertically  as  a  result  of  turbulent  diffusion.  The  horizontal 
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diffusion  is  based  upon  the  commonly  assumed  four-thirds  power  law.  There¬ 
fore,  the  diffusion  coefficient,  (up  to  a  maximum  value  of  100  ft^/s)  is 
given  as 


=  AjL 


4/3 


(125) 


where  is  an  input  dissipation  parameter  and  L  is  set  equal  to  four  standard 
deviations.  As  illustrated  in  Figure  2.4  of  the  report  by  Brandsma  and  Divoky 
(1976),  a  value  of  100  ft^/sec  for  the  horizontal  diffusion  coefficient  cor¬ 
responds  to  a  length  scale  of  lO^-lO'*  ft.  With  the  computational  grid  cell 
typically  being  on  the  order  of  100-500  ft,  a  length  scale  greater  than  1,000  ft 
would  normally  be  associated  with  mean  flow  rather  than  turbulence.  Thus, 
restricting  the  diffusion  coefficient  to  be  less  than  100  ft^/sec  is  reasonable. 

Horizontal  growth  is  achieved  by  employing  the  Fickian  expression 

a„  -  (2K,jr  (1“) 


where 

Ojj  2  =  a  standard  deviation 

t  =  time  since  formation  of  the  cloud 

From  Equation  126. 


(127) 


and  thus, 


a 


=  a 


(128) 


where 


=  a,, ,  at  the  current  time  step 

•^»^new 

a.., ,  ,  ,  =  a., ,  at  the  previous  time  step 
In  a  similar  manner,  the  vertical  growth  is  written  as 
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(129) 


where  is  a  function  of  the  stratification  (including  the  effect  of  the  total 
suspended  sediment)  of  the  water  column.  The  maximum  value  of  Ky  is  input 
as  a  model  coefficient  and  occurs  when  the  water  density  is  uniform. 

If  long-term  output  is  desired  at  the  end  of  a  particular  time-step,  the  con¬ 
centration  of  each  solid  type  is  given  at  each  grid  point  by  summing  the  con¬ 
tribution  from  individual  clouds  to  yield 


C,  -  (2t)-“  Y, 


iri: 


“  ^  ^ 


exp 


-1/2 


(130) 


where  N  is  the  number  of  small  clouds  of  a  particular  solid  type  and  y  (the 
vertical  position  at  which  output  is  desired)  is  specified  through  input  data. 
(x,z)  are  the  horizontal  coordinates  of  a  grid  point.  This  approach  for  the 
transport-diffusion  phase  follows  the  work  of  Brandsma  and  Sauer  (1983). 

The  surface  and  all  solid  boundaries  except  the  bottom  are  handled  by 
assuming  reflection  from  the  boundaries. 

In  addition  to  the  horizontal  advection  and  diffusion  of  material,  settling  of 
the  suspended  solids  also  occurs.  Therefore,  at  each  net  point  the  amount  of 
solid  material  deposited  on  the  bottom  and  a  corresponding  thickness  are  also 
determined.  Since  a  normal  distribution  is  assumed  for  material  in  the  small 
clouds,  deposited  material  is  also  assumed  to  take  such  a  distribution 
horizontally  on  the  bottom.  A  basic  assumption  in  the  model  is  that  once 
material  is  deposited  on  the  bottom,  it  remains  there;  i.e.,  neither  erosion  nor 
bed-load  movement  of  material  is  allowed.  However,  as  previously  discussed, 
deposition  is  prohibited  if  the  computed  bottom  shear  stress  exceeds  a 
specified  critical  shear  stress  for  deposition  for  each  solid  fraction.  This 
allows  for  application  at  dispersive  sites. 

The  discussion  presented  above  for  transport-diffusion  of  solids  also  applies 
to  the  disposed  fluid  with  its  dissolved  constituents.  The  constituents  are 
assumed  to  be  conservative  with  no  further  adsorption  on  or  desorption  from 
the  solids  in  the  water  column  or  those  deposited  on  the  bottom.  Computing 
the  resultant  time-history  of  constituent  concentration  provides  information  on 
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the  dilution  that  can  be  expected  over  a  period  of  time  at  the  disposal  site  and 
enables  the  computation  of  mixing  zones  in  water  column  evaluations. 


Ambient  Environment 

As  illustrated  in  Figure  10,  time-invariant  velocity  profiles  that  allow  for 
flow  reversal  can  be  prescribed.  These  profiles  are  applied  at  each  grid  point. 
Another  option  is  to  specify  a  time-invariant,  spatially  varying  depth-averaged 
velocity.  The  ambient  density  profile  at  the  deepest  point  on  the  grid  must 
also  be  prescribed. 


o)  simple  orthogonal  velocity  profiles  for  constant 

DEPTH.  APPLIED  EVERYWHERE  IN  FIELD. 


b)  VERTICALLY  AVERAGED  VELOCITY  PROFILES  FOR  VARIABLE 
DEPTHS  WITH  EQUIVALENT  LOGARITHMIC  PROFILES 
SUPERIMPOSED. 


Figure  10.  Time  invariant  velocity  profile 


Time-Varying  Faii  Veiocities 

If  a  solid  fraction  is  specified  as  being  cohesive,  the  settling  velocity  is 
computed  as  a  function  of  the  suspended  sediment  concentration  of  that  solid 
type.  The  algorithm  used  is 
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(131) 


=  0.000034  j/  C  s  25  mglt. 

■  =  0.0000225  +  1.6  lo'^  if  25  sCs3000  mg/t 
=  0.0069  if  C  >  3000  mg/l 


where 


=  settling  velocity,  fps 
C  =  suspended  sediment  concentration,  mg/{ 

This  algorithm  was  provided  by  Teeter  (personal  communication)  of  WES. 
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3  Model  Simulations  of 
Laboratory  Tests 


Large-scale  laboratory  tests  of  disposals  from  a  model  split-hull  barge  and  a 
multiple  bin  vessel  are  discussed  by  Johnson,  et  al.  (1993).  The  tests  were 
conducted  in  water  depths  up  to  6  ft  with  the  maximum  horizontal  dimensions 
of  the  test  facility  being  32  ft  by  41  ft.  Both  stationary  and  moving  disposal 
operations  were  simulated  with  disposal  materials  ranging  from  essentially  pure 
clay  to  fine  coal.  Data  collected  consisted  of  bottom  deposition  areas  and 
depths,  suspended  sediment  samples,  and  video  taping. 

Data  from  these  tests  helped  to  guide  the  model  developments  described  in 
Part  2  of  this  report  and  have  also  been  used  to  assess  the  ability  of  STFATE 
to  accurately  simulate  disposal  of  dredged  material  in  open  water.  Results 
from  seven  applications  are  presented  below.  Basic  information  concerning 
the  simulated  tests  are  given  in  Table  1.  For  consistency,  the  test  numbers 
correspond  to  those  reported  in  Johnson,  et  al.  (1993). 


Table  1 

Characteristics  of  Disposal  Tests 

Test  No. 

Type 

Disposal 

Water  Depth 
ft 

Material 

Bulk  Density 
gm/cc 

3 

Barge 

4.5 

Sand 

1.85 

5C 

Barge 

6.0 

Coal 

1.15 

20 

Barge 

6.0 

Clay 

1.13 

22 

Barge 

2.6 

Silt 

1.06 

24 

Barge 

4.0 

Silt 

1.14 

1H 

Hopper 

2.0 

Silt 

1.155 

5H 

Hopper 

4.0 

Silt 

1.18 

As  previously  discussed,  several  model  coefficients  are  required.  Default 
values  were  employed  for  most.  However,  larger  values  were  required  for  the 
descent  and  collapse  entrainment  coefficients.  In  addition,  it  was  found  that 
better  results  were  obtained  if  none  of  the  kinetic  energy  was  expended  to  set 
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the  ambient  fluid  in  motion.  Values  for  all  coefficients  are  given  in  Table  2. 
The  same  values  were  used  in  all  simulations 


Table  2 

Model  Coefficients  used  in  Simulations 

Coefficient 

Description 

Value 

«0 

Descent  Entrainment 

0.60 

p 

Settling 

0.0 

Apparent  Mass 

1.0 

Cd 

Descent  Drag 

1.0 

^DRAG 

Collapse  Drag 

0.5 

^FRIC 

Ruid  Friction 

0.004 

«c 

Collapse  Entrainment 

0.03 

FRICTN 

Bottom  Friction 

0.01 

X 

Horizontal  Diffusion  Parameter 

0.00001 

AKYO 

Unstratified  Vertical  Diffusion  Coeff 

NA 

CSTRIP 

Stripping 

0.003 

An  extremely  low  value  was  used  for  the  horizontal  diffusion  parameter,  X, 
because  the  disposal  environment  was  quiescent.  For  the  same  reason,  the 
value  of  the  vertical  diffusion  coefficient,  AKYO,  is  immaterial.  With  no 
ambient  velocity,  the  laminar  value  for  vertical  diffusion  is  used. 


Test  No.  3 

This  test  was  the  disposal  of  sand  from  a  stationary  1:50  scale  barge  that 
opened  instantaneously.  The  volume  of  material  released  was  0.9  cu  ft, 
representing  about  4000  cu  yd  prototype.  The  disposal  material  had  a  grain 
size  ranging  from  0.15  to  0.80  mm.  As  shown  in  Table  1,  the  water  depth  was 
4.5  ft 

Basic  data  from  the  disposal  test  compared  with  model  results  are  the 
average  descent  velocity  of  the  descenting  jet  or  cloud,  average  rate  of  spread¬ 
ing  of  the  bottom  collapsing  cloud,  depositional  area,  and  suspended  sediment 
concentrations  determined  from  bottle  samples  collected  at  several  depths  at 
various  locations  near  the  barge  at  different  times.  Bottles  placed  on  a  pole 
were  opened  after  the  bottom  surge  front  had  moved  1  ft  past  the  pole.  The 
locations  of  the  poles  containing  the  sediment  bottles  are  shown  in  Figure  11. 
Additional  details  concerning  all  tests  and  the  testing  program  are  provided  by 
Johnson,  et  al.  (1993). 
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Figure  1 1 .  Location  of  sediment  bottles  poles  for  Test  No.  3 


Data  on  the  convective  descent  and  collapse  phases  were  not  available  for 
Test  No.  3.  However,  suspended  sediment  data  were  available  for  compari¬ 
son  with  model  results.  These  are  presented  in  Tables  3-4.  The  suspended 
sediment  concentrations  computed  by  STFATE  at  various  depths  at  increments 
of  1  ft  from  the  point  of  disposal  are  listed.  Measured  concentrations  are 
shown  in  parenthesis  and  were  collected  at  approximately  the  locations 
indicated  on  the  table. 

The  actual  time  at  which  data  were  collected  at  Pole  B  was  not  recorded, 
however,  based  upon  timings  from  other  tests  it  was  probably  in  the  10-15  sec 
range.  Thus,  comparisons  are  shown  at  12  and  15  sec.  Considering  the  fact 
that  the  manner  in  which  the  material  leaves  the  disposal  vessel  and  the  timing 
of  pulling  the  stoppers  from  the  bottles  along  with  the  location  of  the  bottles 
are  crucial  to  the  amount  of  sediment  collected  in  a  bottle,  the  comparison  is 
quite  good.  An  interesting  observation  is  that  in  some  cases  the  measured 
concentrations  were  greater  at  6  in.  off  the  bottom  than  at  3  in.  from  the  bot¬ 
tom.  Table  4  illustrates  that  this  behavior  was  also  computed  by  STFATE. 

An  inspection  of  Tables  3-4  shows  that  the  computed  concentration  field 
changed  significantly  from  12  to  15  seconds.  This  is  because  the  bottom 
surge  is  still  spreading  and  also  because  the  material  stripped  during  descent 
and  from  the  top  of  the  bottom  surge  has  settled  0.3  ft  (assuming  a  fall 
velocity  of  0.1  ^s). 
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Table  3 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/()  at  Pole  B  After  12  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Disposal  (ft) 

0 

1 

2 

3 

4 

5 

6 

■ 

1.5 

3480 

5100 

31 

0 

0 

0 

0 

0 

(1 

1 

N _ 

710 

1640 

270 

3 

0 

0 

0 

I 

(352) 

hz _ 

8700 

18300 

20100 

11900 

4000 

830 

110 

(17074) 

1 

7100 

13400 

17100 

15000 

9000 

3800 

1100 

21^1 

I 

(13907) 

4.25 

11000 

18100 

23800 

25200 

22000 

16700 

11900 

8400 

Table  4 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/C)  at  Pole  B  After  15  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Disposal  (ft) 

0 

1 

2 

3 

4 

5 

6 

7 

1.5 

160 

240 

2 

0 

0 

0 

0 

0 

rznzm: 

(1 

ni 

Li:! _ 

690 

1460 

110 

1 

0 

0 

0 

0  II 

1 

(352) 

lii _ 

6700 

14800 

14400 

630 

130 

10 

0 

1 

(17074) 

6800 

14100 

16100 

10300 

3800 

770 

90 

1 

(13907) 

4.25 

8400 

9700 

10900 

9600 

6400 

3600 

1900 

1200 
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The  measured  concentrations  are  quite  sensitive  to  the  time  at  which  the  bot¬ 
tled  samples  were  taken  since  the  bottom  surge  expands  similar  to  a  circular 
torus.  Samples  taken  in  the  advancing  head  of  the  surge  (see  Figure  12)  reveal 
larger  concentrations  at  greater  heights  above  the  bottom  than  those  taken 
behind  the  head.  Therefore,  it  is  more  meaningful  for  model  and  measured 
results  to  be  compared  in  an  average  sense  over  the  entire  water  column  rather 
than  focusing  on  comparisons  at  a  particular  point  at  a  particular  time.  For 
this  reason,  the  measured  concentrations  for  all  tests  are  presented  between 
concentrations  computed  at  adjacent  depths  and  distances.  The  measured  con¬ 
centrations  should  be  considered  with  respect  to  the  four  computed  values 
which  surrounds  them  on  the  tables. 


Figure  12.  Bottom  surge 


Test  No.  5C 


This  test  was  the  disposal  of  crushed  coal  with  a  density  of  1.3  gm/cc  from 
the  model  split-hull  barge  in  a  water  depth  of  6  ft.  The  average  descent  speed 
of  the  disposal  cloud  was  computed  to  be  0.57  fps.  From  video  tapes  of  the 
disposal,  the  actual  descent  speed  was  0.43  fps.  With  an  entrainment  coeffi¬ 
cient  of  0.60,  the  disposal  material  was  diluted  by  a  factor  of  65  during 
descent.  An  entrainment  coefficient  of  0.60  was  estimated  from  acoustical 
data  collected  during  a  disposal  operation  near  Mobile,  FL  (Kraus  (1991)). 
Based  upon  the  observed  size  of  the  cloud  at  the  moment  of  bottom  impact,  a 
dilution  factor  of  65  is  approximately  correct. 

The  stripping  algorithm  in  STFATE  resulted  in  about  10  percent  of  the 
solids  being  stripped  during  the  descent  through  6  ft  of  water.  Although  no 
measurements  were  taken  that  can  confirm  this,  it  appears  to  be  in  line  with 
published  figures  of  3-5  percent  being  stripped  at  disposal  sites  100-200  ft 
deep.  The  6  ft  water  depth  scales  to  a  prototype  depth  of  300  ft.  At  deeper 
sites,  one  would  expect  a  larger  percentage  of  the  solids  to  be  stripped.  In 


Chapter  3  Model  Simulations  of  Laboratory  Tests 


47 


addition,  the  initial  concentration  of  coal  was  taken  to  be  50  percent,  resulting 
in  a  larger  stripped  quantity. 

The  average  speed  of  the  front  of  the  bottom  surge  over  the  first  3-4  ft 
from  the  point  of  disposal  was  computed  by  STFATE  to  be  0.26  fps.  Results 
from  the  video  data  show  an  advancing  speed  of  0.28  fps. 

Suspended  sediment  comparisons  at  the  locations  shown  in  Figure  13  are 
presented  in  Tables  5-6 
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Figure  13.  Location  of  sediment  botties  poles  for  Test  No.  5C 


Test  No.  20 


This  test  was  the  only  disposal  conducted  with  clay  material.  The  disposal 
occurred  from  the  split-hull  barge  in  6  ft  of  water.  The  averaged  descent 
speed  was  computed  to  be  0.55  fps,  whereas,  the  measured  speed  was  0.67  fps. 
Approximately  2  percent  of  the  solids  are  computed  to  have  been  stripped 
from  the  descenting  cloud.  This  is  significantly  less  than  experienced  by  the 
crushed  coal  disposal  in  Test  No.  5C  due  to  the  cloud  containing  less  sedi¬ 
ment  Stripping  of  material  during  descent  is  related  to  the  concentration  of 
sediment  in  the  descending  cloud. 
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Tabic  5 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/f)  at  Pole  B  After  17  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Disposal  (ft) 

0 

2 

4 

6 

8 

10 

2.5 

425 

360 

0 

0 

0 

0 

1 

(20) 

IU= _ 

3090 

7400 

2060 

70 

0 

1 

(5720) 

1 

in 

500 

1100 

1100 

500 

100 

1 

(18319) 

1 

7800 

17100 

17100 

7800 

1600 

140  1 

1 

(14137) 

1 

1 

5500 

12100 

12000 

5400 

1100 

100 

Table  6 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/f)  at  Pole  A  After  30  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Disposal  (ft) 

0 

2 

4 

6 

8 

10 

12 

14 

1 

(25) 

hi _ 

500 

230 

0 

0 

0 

0 

0 

II 

(216) 

1 

hi _ 

2510 

5790 

1020 

15 

0 

0 

0 

0  1 

1 

(2703) 

4.5 

660 

1580 

800 

90 

2 

0 

0 

0 

(4615) 

II 

5.5 

370 

710 

900 

760 

430 

160 

40 

0  II 

(2347) 

II 

5.75 

1220 

2330 

2970 

2520 

1430 

540 

140 

25 
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The  average  speed  of  the  bottom  surge  over  the  first  couple  of  feet  was 
computed  to  be  0.38  fps,  whereas,  the  speed  determined  from  video  tapes  was 
0.57  fps.  The  reason  for  the  difference  is  not  known.  Model  and  measured 
bottom  surge  speeds  compared  well  in  most  tests.  The  difference  could  be  due 
to  reduced  frictional  effects  at  the  bottom. 


Suspended  sediment  concentrations  at  the  locations  shown  in  Figure  14  are 
presented  in  Tables  7-8. 


- 

4 

ZD  - 

pS- 

^  - 

to 

CO  - 

to 

ffl 

■ 

■ 

i 

■ 

Cfc:r 

X 

cz 

a " 

>-• 

a.  - 

O  " 

^  “ 

21  - 
11 

- r 

1. 

- 1 - 1 

M  I 

X 

- r 

5  I 

RXI5 

r 

2 

% 

I 

3  2. 

1 

2  2 

Figure  14.  Location  of  sediment  bottles  poles  for  Test  No.  20 


Model  results  and  measured  concentrations  compare  well,  although,  the  model 
computes  concentrations  slightly  smaller  than  those  measured  in  the  upper 
water  column. 

Bottom  deposition  was  not  measurable.  However,  the  general  observation 
that  the  material  covered  the  entire  test  area  agrees  with  the  model  results. 
Computed  depositions  greater  than  0.05  x  10'^  ft  covered  an  area  with  a 
diameter  of  32  ft.  The  maximum  deposition  thickness  was  computed  to  be 
3.3  X  10'^  ft. 


Chapter  3  Model  Simulations  of  Laboratory  Tests 


better  agreement  than  those  in  shallower  water  which  implies  that  the 
algorithm  for  computing  the  injection  speed  should  be  investigated. 

The  average  bottom  surge  speed  over  the  first  couple  of  feet  was  computed 
to  be  0.46  fps.  Video  tape  data  shows  an  actual  speed  of  about  0.42  fps. 

Computed  suspended  sediment  concentrations  at  the  locations  shown  in 
Figure  15  along  with  the  measured  concentration  placed  in  parenthesis  are 
presented  in  Tables  9-11. 


Figure  15.  Location  of  sediment  bottles  poles  for  Test  No.  22 


An  inspection  of  these  results  reveals  that  the  model  underpredicted  concentra¬ 
tions  in  the  upper  part  of  the  water  column  as  the  surge  expanded.  It  should 
be  realized  that  the  exact  shape  of  the  bottom  surge  is  not  simulated. 

STFATE  assumes  that  the  surge  is  a  collapsing  half  ellipsoid,  whereas,  the 
actual  surge  is  a  circular  torus.  The  effect  of  the  torus  shape  on  water  column 
concentrations  is  approximately  accounted  for  in  STFATE  by  employing  an 
offset  Gaussian  distribution  for  the  clouds  stripped  away  at  the  top  of  the 
surge.  As  illustrated  in  Tables  9-10,  this  results  in  better  predictions  of 
suspended  sediment  concentrations  in  the  upper  water  column  closer  to  the 
discharge  point. 

The  deposited  material  covered  an  area  of  about  200  fl^  with  an  elliptical 
shape.  The  computed  area  of  deposition  with  deposits  of  greater  depth  than 
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Table  9 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/f)  at  Poie  B  After  12  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Disposal  (ft) 

0 

1 

2 

3 

4 

5 

6 

0.5 

405 

275 

0 

0 

0 

0 

0 

(429) 

1.5 

4170 

8650 

520 

0 

0 

0 

(13C») 

2.0 

3600 

8730 

4410 

460 

10 

0 

j _ 

(3306) 

2.25 

2310 

5100 

— 

5040 

2250 

470 

50 

5 

Tabie  10 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/C)  at  Poie  A  After  23  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Discharge  (ft) 

0 

1 

2 

3 

4 

5 

0.5 

210 

170 

0 

0 

0 

(288)  1 

1.5 

3380 

6880 

360 

360 

0 

r 

(615)  1 

3430 

8280 

2520 

2520 

0 

0 

(388)  1 

2.25 

2530 

5760 

5070 

1720 

230 

10 

0.03  X  10'^  ft  was  about  150  ft^,  with  the  greatest  deposition  depth  being  1.9  x 
10'^  ft. 


Test  No.  24 

This  test  was  the  disposal  of  a  silt  slurry  with  a  density  of  1.14  gm/cc  in 
4.0  ft  of  water  from  the  model  split-hull  barge.  The  average  descent  speed 
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Table  1 1 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/f )  at  Pole  C  After  37  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Discharge  (ft) 

0 

1 

2 

3 

4 

5 

0.5 

45 

50 

0 

0 

0 

0 

(149) 

1.5 

195 

370 

20 

0 

0 

0 

(423) 

2.0 

3520 

8440 

2260 

70 

— 

0 

0 

(609) 

2.25 

2830 

6760 

3960 

730 

75 

5 

was  computed  to  be  0.80  fps,  whereas,  the  measured  speed  was  0.73  fps. 

The  dilution  factor  was  computed  to  be  about  20  with  approximately  1.5  per¬ 
cent  of  the  solids  stripped  during  descent. 

The  average  speed  of  the  front  of  the  bottom  surge  over  the  first  couple  of 
feet  was  computed  to  be  0.42  Ips  with  the  measured  speed  being  0.28  fps.  As 
previously  noted,  material  can  leave  the  disposal  vessel  in  a  nonuniform 
fashion  rather  than  as  a  single  cloud  formed  instantaneously.  This  can  have 
an  impact  on  the  speed  of  the  surge  during  its  early  development. 

A  comparison  of  computed  and  measured  suspended  sediment  concentra¬ 
tions  at  the  locations  shown  in  Figure  16  is  given  in  Tables  12-14. 

Unlike  some  of  the  previous  comparisons,  these  results  show  better  agree¬ 
ment  away  from  the  disposal  point.  Again,  the  exact  reason  is  not  known  but 
as  previously  discussed,  the  time  of  the  release  of  the  bottle  corks  and  the 
manner  in  which  material  is  released  play  a  major  role  in  the  amount  of 
material  captured  by  an  individual  bottle. 

The  observed  area  of  deposition  was  about  225  fF  and  was  elliptical  in 
shape.  The  computed  area  of  deposition  within  which  the  depth  exceeded 
0.21  X  10'^  ft  was  about  200  fF  with  the  maximum  depth  being  1.60  x  10'^  ft. 
Depths  less  than  0.21  x  10'^  ft  were  outside  the  computational  grid. 


Test  No.  1H 


This  test  was  the  disposal  of  a  silt  slurry  from  a  1:50  scale  model  hopper 
vessel  containing  6  bins.  The  bulk  density  was  1.155  gm/cc  and  the  total 
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Figure  16.  Location  of  sediment  bottles  poles  for  Test  No.  24 


volume  was  1.68  ft^.  Disposal  occurred  in  2.0  ft  of  water  with  pairs  of  bins 
opened  in  sequence.  With  the  disposal  occurring  in  such  shallow  water, 
impact  with  the  bottom  occurred  almost  instantaneously.  The  average  speed  of 
the  surge  front  over  the  first  2-3  ft  was  0.44  fps  with  the  average  measured 
speed  being  0.42  fps. 

Comparisons  of  computed  and  measured  suspended  sediment  concentrations 
at  the  locations  shown  in  Figure  17  are  given  in  Tables  15-16.  In  these  tables, 
four  rows  of  computed  concentrations  are  presented  at  each  depth  with  the  first 
row  lying  along  the  center  line  of  the  vessel. 


Test  No.  5H 

This  test  was  similar  to  Test  IH.  However,  the  disposal  was  in  4.0  ft  of 
water  and  the  bulk  density  of  the  silt  slurry  was  1.18  gm/cc.  The  average 
computed  descent  speed  was  1.10  fps,  whereas,  the  measured  speed  was 
1.48  fps.  The  average  speed  of  the  bottom  surge  over  the  first  2-3  ft  was  com¬ 
puted  to  be  0.47  fps  with  the  average  measured  speed  being  0.45  fps. 

Comparisons  of  computed  and  measured  suspended  sediment  concentrations 
at  locations  shown  in  Figure  18  are  presented  in  Tables  17-18.  As  in 
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Table  12 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/f)  at  Pole  B  After  12  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Discharge  (ft) 

0 

1 

2 

3 

4 

5 

6 

■ 

8 

1.0 

640 

870 

5 

0 

0 

0 

0 

2.0 

490 

1160 

490 

55 

0 

0 

0 

■1 

Hi 

(562) 

3.0 

1330 

2870 

2990 

1550 

410 

55 

5 

■9 

0 

(1544) 

1440 

2580 

3410 

3310 

2370 

1250 

490 

150 

30 

(1932) 

_ ^1 

3.75 

2530 

4280 

5740 

6110 

5140 

3430 

1820 

760 

250 

Table  13 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/f)  at  Pole  A  After  21  Seconds 

Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Discharge  (ft) 

0 

1 

2 

3 

4 

5 

6 

■ 

8 

1.0 

355 

500 

0 

0 

0 

0 

0 

^9 

- 

(129) 

2.0 

165 

380 

50 

0 

0 

0 

0 

0 

(918) 

3.0 

1930 

4000 

4550 

2890 

1030 

200 

20 

m 

0 

(2353) 

3.5 

1730 

3190 

4160 

3820 

2490 

1150 

370 

90 

0 

(1900) 

3.75 

680 

1020 

1350 

1590 

1660 

1530 

1260 

920 

590  1 
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Table  15 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/()  at  Pole  B  After  15  Seconds 


Horizontal  Position  Relative  to  Point  of  Discharge  (ft) 

Depth 

(ft) 

0 

1 

2 

3 

4 

5 

6 

1.0 

14900 

21500 

200 

0 

0 

0 

0 

21500 

5800 

35 

0 

0 

0 

0 

200 

35 

0 

0 

0 

0 

0 

(529) 

0 

0 

0 

0 

0 

0 

0 

1.5 

73000 

104000 

4700 

80 

0 

0 

0 

104000 

33000 

2000 

30 

0 

0 

0 

4700 

2000 

170 

0 

0 

0 

0 

(421) 

80 

30 

0 

0 

0 

0 

0 

Depth 

(ft) 

0 

1 

2 

3 

4 

5 

6 

1.75 

29100 

58400 

18100 

3300 

210 

5 

0 

58400 

37200 

13200 

4800 

710 

40 

0 

18100 

13200 

4800 

710 

40 

0 

0 

(1806) 

1 

3300 

2300 

710 

90 

5 

0 

0 

Table  16 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/£)  at  Pole  A  After  27  Seconds 

Depth 

(ft) 

Horizontal  Position  Reiative  to  Point  of  Discharge  (ft) 

0 

1 

2 

3 

4 

1.0 

1960 

540 

0 

0 

0 

540 

25 

0 

0 

0 

0 

0 

0 

0 

0 

(1271) 

(Continued) 
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Table  17 

Comparison  of  Computed  and  Measured  Suspended  Sediment 
Concentration  (mg/()  at  Pole  C  After  1 1  Seconds _ 


Depth 

(ft) 

Horizontal  Position  Relative  to  Point  of  Discharge  (ft)  || 

0 

1 

2 

3 

4 

5 

6 

1.0 

1140 

860 

0 

0 

0 

0 

0 

860 

90 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(737) 

1 

0 

0 

0 

0 

0 

0 

0 

2.0 

6500 

15300 

6500 

540 

10 

0 

0 

15300 

13000 

4200 

310 

5 

0 

0 

6500 

4200 

930 

55 

0 

0 

0 

(203) 

540 

310 

55 

0 

0 

0 

3.0 

10800 

25100 

12700 

3000 

370 

30 

0 

25100 

21800 

9500 

2200 

280 

20 

0 

12700 

9500 

4000 

910 

110 

10 

0 

(738) 

3000 

2200 

910 

210 

30 

0 

0 

3.75 

8000 

16900 

17600 

10600 

4400 

1300 

0 

16900 

18600 

16200 

9400 

3800 

1100 

0 

17600 

16200 

11900 

6500 

2500 

700 

0 

(4091) 

_ 

10600 

9400 

6500 

3300 

1300 

340 

0 

Tables  15-16,  several  rows  of  computed  concentrations  are  presented  at  each 
depth.  It  can  be  seen  that  computed  and  measured  results  usually  agree  quite 
well  within  one  grid  point. 

The  observed  bottom  deposition  covered  virtually  the  entire  test  section. 
Model  results  substantiated  this  with  the  maximum  computed  thickness  being 
4.2  X  10’^  ft. 
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4  Summary  and  Conclusions 


Summary 

Numerical  models  are  required  to  assess  water  column  environmental 
impacts  of  dredged  material  disposals  in  open  water.  Such  short-term  fate 
models  also  provide  estimates  of  the  initial  deposition  of  sediments  that  are 
then  used  in  determining  the  long-term  movement  of  material  within  the 
disposal  site.  The  earliest  development  of  a  comprehensive  numerical  model 
for  predicting  the  physical  fate  of  material  disposed  in  the  ocean  was  the  work 
of  Koh  and  Chang  (1973).  However,  even  though  additional  developments  of 
the  Koh  and  Chang  model  continued  over  the  next  several  years,  deficiencies 
still  remained  at  the  initiation  of  the  US  Army  Corps  of  Engineers  Dredging 
Research  Program  (DRP)  in  1988. 

To  provide  guidance  on  removing  these  deficiencies,  large  scale  (1:50) 
laboratory  tests  of  split-hull  barge  and  hopper  dredge  disposals  were  conducted 
under  the  DRP.  Results  from  these  tests  are  presented  by  Johnson,  et  al. 
(1993).  The  numerical  model  that  has  evolved  from  the  DRP  is  called 
STFATE  (Short-Term  FATEl.  In  addition  to  guiding  model  developments, 
results  from  the  laboratory  tests  have  been  employed  to  assess  the  ability  of 
STFATE  to  accurately  simulate  open  water  disposal  operations. 

Efforts  to  remove  model  deficiencies  included  developments  in  the  follow¬ 
ing  areas:  (a)  computation  of  insertion  velocity,  (b)  dynamic  collapse  on  the 
bottom,  (c)  effect  of  bottom  slope,  (d)  disposal  at  dispersive  sites,  (e)  consoli¬ 
dation  and  hopper  dredge  disposal,  (f)  stripping  of  sediment  and  fluid, 

(g)  time-varying  fall  velocities,  and  (h)  distribution  of  bottom  deposits.  These 
developments  have  resulted  in  a  more  realistic  simulation  of  real  disposal 
operations.  For  example,  the  use  of  multiple  convecting  clouds  allows  for  a 
realistic  simulation  of  disposal  from  hopper  dredges  and/or  the  effect  of 
consolidation  of  material  during  transit  to  the  disposal  site. 

Results  from  model  simulations  of  seven  disposal  tests  have  been  presented. 
These  include  five  disposals  from  the  model  split-hull  barge  and  two  from  the 
multiple-bin  disposal  vessel.  Water  depths  ranged  from  2.0  to  6.0  ft  and 
sediment  types  were  sand,  silt,  clay,  and  fine  crushed  coal.  At  a  1  to  50  scale, 
these  tests  represented  disposal  of  about  4,000  cu  yd  from  the  split-hull  barge 
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and  8,000  cu  yd  from  the  hopper  dredge  type  vessel  in  water  depths  of  100- 
300  ft. 


Conclusions 

From  an  inspection  of  the  model  results  compared  with  measured  data  from 
the  disposal  tests,  it  can  be  concluded  that  STFATE  simulates  the  fate  of 
material  disposed  in  open  water  well.  Average  computed  and  measured 
descent  speeds  compared  to  within  25  percent  and  average  bottom  surge 
speeds  compared  to  within  about  10  percent.  Considering  the  idealized  geo¬ 
metrical  shapes  assumed  in  STFATE  and  the  uncertainty  of  material  char¬ 
acteristics  and  the  manner  in  which  material  leaves  the  disposal  vessel,  i.e.  the 
initial  conditions,  these  comparisons  seem  reasonable. 

All  of  the  numerical  model  simulations  allowed  sediment  to  be  stripped. 
Computed  results  for  the  silt  and  clay  disposals  showed  about  2-3  percent  of 
the  solids  being  stripped  during  descent.  These  results  agree  with  published 
estimated  percentages  of  stripped  material. 

When  comparing  computed  and  measured  suspended  sediment  concentra¬ 
tions  collected  by  the  sediment  bottles,  it  should  be  realized  that  small  differ¬ 
ences  in  timing  and  the  location  of  the  bottles  are  extremely  important.  In 
addition,  the  manner  in  which  material  leaves  the  vessel  is  important.  In  some 
tests,  small  "globs"  of  material  left  the  vessel  after  the  main  body  of  the 
descending  jet  had  impacted  the  bottom.  Although  the  concept  of  multiple- 
convecting  clouds  allows  for  the  effect  of  consolidation,  the  behavior  described 
above  is  impossible  to  predict. 

Suspended  sediment  results  from  the  disposal  of  sand  in  4.5  ft  of  water 
from  the  model  barge  compared  well  with  computed  results  near  the  bottom. 
Both  computed  and  measured  concentrations  were  about  17,000  mg/C.  How¬ 
ever,  computed  concentrations  were  too  high  in  the  upper  water  column. 
Perhaps  stripping  of  sand  should  not  have  been  allowed. 

Computed  concentrations  of  crushed  coal  disposed  in  6.0  ft  of  water  com¬ 
pared  well  at  virtually  all  locations.  Near  the  point  of  discharge,  measured 
concentrations  varied  from  about  20  to  18,000  mg/C  over  the  water  column. 
Computed  values  ranged  from  near  zero  to  about  17,000  mg/C.  At  a  location 
farther  from  the  discharge,  measured  concentrations  ranged  from  25  to 
4,600  mg/C,  whereas,  computed  values  ranged  from  zero  near  the  surface  to 
3,000  mg/C  near  the  bottom. 

Computed  and  measured  concentrations  for  the  clay  disposal  in  6.0  ft  of 
water  compared  extremely  well  throughout  the  water  column  away  from  the 
point  of  discharge.  These  concentrations  were  in  the  50-1,000  mg/C  range. 

Two  simulations  of  the  disposal  of  a  silt  mixture  from  the  model  barge  in 
water  depths  of  2.6  and  4.0  ft  of  water  were  conducted.  Comparisons  of 


Chapter  4  Summary  and  Conclusions 


concentrations  were  excellent  close  to  the  discharge  point  for  disposal  in  2.6  ft 
of  water  with  concentrations  in  the  500-3,500  mg/{  range.  However,  away 
from  the  disposal  site,  model  results  were  too  low  in  the  upper  water  column. 
Comparisons  for  the  4.0-ft  disposal  were  opposite  from  those  for  the  disposal 
in  2.6  ft  of  water.  Model  results  showed  concentrations  ranging  from  500- 
6,000  mg/{  near  the  discharge  point  with  measured  data  being  in  the  500- 
2,000  mg/J  range.  At  a  location  slightly  farther  away,  both  computed  and 
measured  concentrations  were  in  the  1,000-2,500  mg/{  range.  These  differ¬ 
ences  can  be  attributed  to  the  circular  torus  shape  of  the  bottom  surge  and  the 
manner  in  which  STFATE  attempts  to  model  it. 

Results  from  the  multi-bin  disposals  showed  that  computed  results  generally 
compared  well  with  measured  values  near  the  bottom,  but  were  too  low  in  the 
upper  water  column.  This  was  true  for  disposal  in  both  2.0  ft  and  4.0  ft  of 
water.  Again,  it  appears  the  shape  of  the  surge  is  responsible  for  the 
differences. 

For  most  of  the  disposal  tests,  a  comparison  of  the  bottom  deposition  could 
only  be  accomplished  in  a  qualitative  sense  since  the  depth  of  deposition  was 
too  small  to  measure.  However,  model  results  generally  agreed  well  in  a 
qualitative  sense  with  observed  areas  of  deposition. 

Results  from  these  simulations  have  substantiated  that  STFATE  can  be  used 
to  accurately  simulate  the  fate  of  material  during  disposal  operations.  Descent 
and  bottom  surge  speeds,  stripping,  rates  of  dilution,  total  depositional  areas 
and  suspended  sediment  concentrations  in  the  bottom  surge  are  all  reasonably 
reproduced.  However,  it  is  recommended  that  further  attention  be  given  to 
more  accurately  represent  the  effect  of  the  torus  shape  of  the  bottom  surge. 

An  example  of  additional  research  that  is  needed  to  make  STFATE  even 
more  useful  for  addressing  environmental  issues  lies  in  the  area  of  uncertainty. 
Given  the  uncertainty  in  specifying  characteristics  of  the  disposal  material  at 
the  moment  of  disposal,  the  manner  in  which  material  leaves  the  disposal 
vessel,  and  ambient  conditions,  an  uncertainty  analysis  should  be  conducted  to 
better  define  bounds  on  the  accuracy  of  predicted  water  column  effects. 
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